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ABSTRACT 

Using a statistical sample of dark matter haloes drawn from a suite of cosmological iV- 
body simulations of the Cold Dark Matter (CDM) model, we quantify the impact of a sim- 
ulated halo's mass accretion and merging history on two commonly used measures of its 
dynamical state, the virial ratio -q and the centre of mass offset Ar. Quantifying this relation- 
ship is important because the degree to which a halo is dynamically equilibrated will influence 
the reliability with which we can measure characteristic equilibrium properties of the struc- 
ture and kinematics of a population of haloes. We begin by verifying that a halo's formation 
redshift Zform correlates with its virial mass Afvir and we show that the fraction of its recently 
accreted mass and the likelihood of it having experienced a recent major merger increases 
with increasing Afvir and decreasing Zform- We then show that both 77 and Ar increase with 
increasing Mvir and decreasing Zform^ which implies that massive recently formed haloes are 
more likely to be dynamically unrelaxed than their less massive and older counterparts. Our 
analysis shows that both 77 and Ar are good indicators of a halo's dynamical state, showing 
strong positive correlations with recent mass accretion and merging activity, but we argue 
that Ar provides a more robust and better defined measure of dynamical state for use in cos- 
mological iV-body simulations at z ~ 0. We find that Ar < 0.04 is sufficient to pick out 
dynamically relaxed haloes at z=0. Finally, we assess our results in the context of previous 
studies, and consider their observational implications. 

Key words: methods: iV-body simulations - galaxies: formation - galaxies: haloes - cos- 
mology: theory - dark matter - large-scale structure of Universe 



1 INTRODUCTION particularly important one when seeking to characterise the struc- 
ture and kinematics of simulated haloes in cosmological simu- 

One of the fundamental assumptions underpinning modem theories lations. Previous studies have shown that dynamically unrelaxed 

of galaxy formation is that gal axies form and evolve in massive viri- haloes tend to have lower central densities (see, for example , 

alised haloes of dark matter t White & Ree sl l 19781 : Iwhite & Fren3 IXormen et al.ll 19971 ; iMaccio et alj|2007l : iRomano-Diaz et alj|2007h 

Il99lh . Characterising the properties of these haloes is an impor- and h i gher velocity dispersions (se e , for example, IXormen et all 

tant problem, both theoretically and observationally, and its study 1 19971 : faetznecker & Burkerll l2006l : iD'Onghia <fe Navarro! l2007h 



has been one of the main objectives of cosmological iV-body sim- than their dynamically relaxed counterparts. This means that a dy- 

ulations over the last two decades. The majority of these Simula- namically unrelaxed halo is likely to have a measurably lower con- 

tions have modeled halo formati on and evolution in a purely Cold centration Cvir and higher spin para meter A than its dynamically 

Dark Matter (CDM) universe (cf. lSpringel et alj2006h, with the fo- relaxed counterpart (see, for example. lGardneJl200ll : lMacci6 et all 

cus primarily on their equilibrium structure (cf.'Diemand & Moord |2007h . and so care must be taken to avoid contaminating halo sam- 

L2OO9). Various studies have revealed that CDM hal oes in dy- pies with dynamically unrela xed systems whe r i measuring, for ex 



nami cal equilibrium are triaxial structures (e.g. iBailin & Steinmetzl ample, spin distributions (e g.lSett et alJl20()7l:lMacci6 et alj|2007l : 



20051) supporte d by velocity dispersion rather than rotation (e.g. iD'Onghia & Navarrol2007l:lKnebe & Powej2008l) and the correla- 



Bett et al.l |2007b, with mass pro files that are divergen t down to the tion of halo mass and concentration Cyi r - M,,w (e.g.lMaccio et al.1 

smallest resolvable radius (e.g. 'Diemand et alj|2008|; IStadel et all |2007l : lNeto et alj2007l : lGao et alj2008l : |Prada et alj201lh . 

2009; Navarro et^ 2010) and an abundance of substructure (e.g. yet haloes do not exist in isolation, and the degree to which 

Diemand et aH2007| ;| Sprmgel et al. || 2008| ;| Gao et al. || 201 1 ] ) . ^^^^ dynamically relaxed or unrelaxed bears the imprint of 



The qualification that a halo is in dynamical equilibrium is a both their environment and their recent mass assembly and merg- 
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ing history. As previous studies have shown, dynamically unrelaxed 
haloes t end to have suffered one or more recent signific ant merg- 
ers (e.g. lTormen et al.lll997l : lHetznecker & Burker For this 
reason, it is common practice to use dynamical state and recent 
merging history interchangeably, with the understanding implicit 
that unrelaxed haloes are ones that have suffered one or more re- 
cent major mergers. 

However, it is important to establish this practice on a more 
quantitative footing and to assess how well a halo's dynamical state 
and its recent mass assembly history correlate. This is because of 
the need to identify robustly haloes that are in dynamical equilib- 
rium - or indeed disequilibrium - in cosmological simulation^ 
The goal of this paper is to quantify this relationship using a sta- 
tistical sample of haloes drawn from cosmological A^-body simu- 
lations of the CDM model. The CDM model is the ideal testbed for 
this study because of the fundamental role merging plays in halo 
mass assembly (e.g. Maulbetsch et al. 2007; Fakhouri & Ma 2008,; 
iMcBride et alll2009l ; iFakhouri & M^|2010|; iFakhouri et al.ll2010h . 
and because we expect massive haloes, which on average form later 
than their less massive counterparts, to have more violent recent 
merging histories. 

Such an undertaking has practical implications. For example, 
if we want to robustly characterise the predicted variation of, say, 
concentration Cvir with virial mass Afvir on galaxy group and clus- 
ter mass scales (Mvir ^ IO^^Mq), then it is essential that we can 
identify relaxed systems in a robust fashion. Should we use mass 
assembly histories directly and select only haloes that have quies- 
cent recent merging histories, or are commonly used measures that 
estimate dynamical state based on material within the halo's virial 
radius rvir adequate? This is particularly important for comparison 
with observations that provide crucial tests of the theory, such as 
the analysis the A/vir — Cvir relation for groups and clusters dra wn 
from the Sloan Digital Sky Survey bv lMandelbaum et al.] j2008h . 

In this paper, we examine how a halo's mass assembly history 
and dynamical state varies with its virial mass A/vir and its for- 
mation redshift, and adopt simple measures to characterise a halo's 
recent mass assembly and merging history - namely, the fraction of 
mass assembled (AM/M); the rate of change of mass with redshift 
1/MdM/dz; and the most significant merger 5max. We compare 
these with two measures of the halo's dynamical state - the virial 
ratio 

r,^2T/\W\, (1) 

where T and W are t he kinetic and gravitational potential energies 
of ha lo material (cf. ICole & Lacev. .1996, ; .Hetznecker & BurkertI 
l2006l) . and the centre-of-mass offset 

Ar = |fccn — rem] /rvir, (2) 

where fccn and fcm are the centres of density and mass of halo 
material and rvir is t he hal o's virial radius (cf. ICrone et alj[r996l ; 
iThomas et al. I 1 19981 |2001[) . Previous studies have shown that 
both Tj and Ar increase in the aftermath of a major merger (e.g. 
iHetznecker & BurkerL2006. ; .Poole et al.,,2006,) , and we will clarify 

^ Our focus is fixed firnily on haloes in cosmological simulations, but 
we note that the relationship between dynamical state and recent mass 
assembly history is equally important observationally. Here, for exam- 
ple, estimates of the dynamical masses of galaxy clusters require as- 
sume a p opulation of dynamical trac ers that are in dynamical equilib- 
rium (e.g. IPiffaretti & Valdaminill20o3) . while reconstructions of a galaxy 
cluster's recent mer ging history look for signatures of disequilibrium (e.g. 
ICassano et alj20l"ol) . See §[6]for further discussion. 



precisely how they relate to a halo's mass assembly and merging 
activity in g eneral. We note that ou r work develops earlier ideas 
presented in iKnebe & Powe3 i l2008l) . in which we investigated the 
relationship between halo mass A/vir and spin A, and it comple- 
ments that of lOavis et al. Ilmil), who address related but distinct 
issues in their critique of the application of the virial theorem (cf. 
^4.11 to simulated high redshift dark matter haloes. 

The layout of the paper is as follows. In ^ we describe our 
approach to making initial conditions; finding and analysing dark 
matter haloes in evolved outputs; constructing merger trees of our 
dark matter haloes; and our criteria for defining our halo sample. 
In !j3] we examine the relationship between a halo's virial mass 
A/vir, its formation time Zform and measures of its mass accretion 
and merging history. In ^ we present commonly used measures 
for assessing the dynamical state of a dark matter halo - the virial 
ratio -q — 2T/\W\ (cf. ^4Ai and the centre-of-mass offset Ar = 
|rcen — rem I /rvir (cf. i|4.2| l - and investigate how these measures 
correlate with A/vir and Zform- In !j5]we combine the insights from 
the previous two sections and show how a halo's dynamical state 
depends on its recent mass accretion and merging history. Finally, 
we summarise our results in ij6]and comment on the implications 
of our findings for both observational studies and galaxy formation 
modeling. 



2 METHODS 

2.1 The Simulations 

We have run a series of 256"^ particle cosmological A^-body sim- 
ulations following the formation and evolution of structure in the 
CDM model. We use a sequence of boxes of side Lbox varying be- 
tween 20/i^^Mpc and 500/i^^Mpc from Zstart=100 to Zflnish=0. 
In each case we assume a flat cosmology with a dark energy term, 
with cosmological parameters fio ~ 0.7, JIa = 0.3, h — 0.7, and 
a normalisation as — 0.9 at 2: = 0. Various properties of these 
simulations are summarised in table[T] 

Initial conditions were generated using a standard procedure 
that can be summarised as follows; 

(i) Generate the CDM transfer function for the appropriate cos- 
mological param eters (Qp , Qa , and h) u sing the Boltzmann 
code CMBFAST jSeliak & Zaidarriagal ri996h . This is convolved 
with the primordial power spectrum P{k) oc k", n = 1, to obtain 
the unnormalised power spectrum, which is normalised by requir- 
ing that the linear mass a{R) equal as on a scale of 8/i~^Mpc at 

2=0. 

(ii) Create a statistical realisation of a Gaussian random field of 
density perturbations in Fourier space, whose variance is given by 
P{k), where k = ^^k'^ + ky + k^ and whose mean is zero. 

(iii) Take the inverse transform of the density field and compute 
positions and velocities using the Zel'dovich approximation. 

(iv) Impose these positions and velocities on an initial uniform 
particle distribution such as a grid or "glass". 

Note that throughout our we use a " glass"-like co nfiguration as our 
initial uniform particle distribution ( IWhite|[r996h . 

A ll simulations w ere run using the parallel TreePM code GAD- 
GET! (ISpringelll2005l) with constant comoving gravitational soft- 
ening £ and individual and adaptive timesteps for each particle. 
At = r\\fTfa, where a is the magnitude of a particle's gravita- 
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Table 1. Properties of the Simulations. Each of the simulations contains 
256^ particles. In addition, L^ox is the comoving box length in units of 
/i~^Mpc; Njun is the number of runs in the series; nipart is the particle 
mass in units of /i^^Mq; e is the force softening in comoving units of 
/i~^kpc; and Afcut is the halo mass corresponding to A^cut=600 particles, 
in units of h~^MQ- 



Run 


-^'box 


Nrun ^part 


e 


A/cut 


ACDMX20 


20 


5 3.97 X lO'^ 


1.5 


2.38 X 10^" 


ACDMX50 


50 


1 6.20 X 10* 


3.9 


3.72 X 10" 


ACDMX70 


70 


1 1.70 X 10^ 


5.5 


1012 


ACDMXIOO 


100 


1 4.96 X 10^ 


7.8 


2.97 X 10^2 


ACDMX200 


200 


1 3.97 X 10^° 


15.6 


2.39 X 10^^ 


ACDMX500 


500 


1 6.20 X 10" 


39.1 


3.72 X lO^* 



tional acceleration and rj = 0.05 determines the accuracy of the 
time integration. 



2.2 Halo Identification and Merger Trees 

Halo Catalogues Groups were identified using the MPI-enabled 
version of AHF, otherwise known as AMIGA'S Halo Finde{3 
jKnoUmann & Knebd l2009h. AHF is a modification of MHF 



(MLAPM's Halo Finder; see'Gill et aL''2004') that locates groups as 
peaks in an adaptively smoothed density field using a hierarchy of 
grids and a refinement criterion that is comparable to the force res- 
olution of the simulation. Local potential minima are calculated for 
each of these peaks and the set of particles that are gravitationally 
bound to the peaks are identified as groups that form our halo cata- 
logue. 

For each halo in the catalogue we determine its centre-of- 
density fc cn (using the itera tive "shrinking spheres" method de- 
scribed in IPower et al.ll2003h and identify this as the halo centre. 
From this, we calculate the halo's virial radius rvir, which we 
define as the radius at which the mean interior density is Avir 
times the critical density of the Universe at that redshift, Pc{z) = 
3H'^{z)/8tvG, where H{z) and G are the Hubble parameter at 
z and the gravitational constant respectively. The corresponding 
virial mass A/vir is 



Afvir 



47r • 3 

-^ZAvirPc^vir- 



(3) 



We adopt a cosmology- and redshift-dependent overdensity crite- 
rion, which for a ACDM cosmology with Qo = 0.3 and Qa = 0.7 
gives Avir 97 at z=0 (c.f. iEke et al...l99a) . 

Merger Trees Halo merger trees are constructed by linking halo 
particles at consecutive output times; 

• For each pair of group catalogues constructed at consecutive 
output times ti and t2 > ti, the "ancestors" of "descendent" 
groups are identified. For each descendent identified in the cata- 
logue at the later time t2, we sweep over its associated particles 
and locate every ancestor at the earlier time ti that contains a sub- 
set of these particles. A record of all ancestors at ti that contain 
particles associated with the descendent at t2 is maintained. 



AHF may be downloaded from http : / /popia . ft . uam . es /AMIGA 



• The ancestor at time ti that contains in excess of /prog of these 
particles and also contains the most bound particle of the descen- 
dent at t2 is deemed the main progenitor. Typically ./prog ~ 0.5, 
i.e. the main progenitor contains in excess of half the final mass. 

Each group is then treated as a node in a tree structure, which 
can be traversed either forwards, allowing one to identify a halo 
at some early time and follow it forward through the merging hi- 
erarchy, or backwards, allowing one to identify a halo and all its 
progenitors at earlier times. 



2.3 Defining the Halo Sample 

A degree of care must be taken when choosing which haloes to in- 
clude in our sample, to ensure that our results are not affected by the 
finite resolution of our simulations. One of the key calculations in 
this study is of a halo's virial ratio rj = 2T/\W\ (see ij |4.1| l, where 
T and W are the kinetic and gravitational potential energies of ma- 
terial within Tvir- The gravitational potential energy is particularly 
sensitive to resolution; if a halo is resolved with too few particles, 
its internal structure will not be recovered sufficiently accurately 
and the magnitude of W will be underestimated. 

We estimate how many particles are needed to recover W ro- 
bustly from a TV-body simulation in Figure [T] Here we generate 
Monte Carlo A'^-body realisations of a halo whose spherica lly aver- 
aged mass profile is described by the lNavarro et alj ( Il997b profile, 



Pc 



cx (1 -|- cx) 



2 ' 



(4) 



here x = r/rvir is the radius r normalised to rvir, c is the concen- 
tration parameter and Sc is the characteristic density, 

c ^ Avir 

3 ln(l + c) - c/(l + c) ■ ^ ' 

The resulting gravitational potential energy is given by 

"c (2 + c) ln(l + c) 



W = 



-16n^ G Pc 5c 



2(l + c)2 (1 + c) 



(6) 



In a TV-body simulation or realisation, we calculate W by randomly 
sampling particles within rvir and reseating; this gives 



W - 



■ATv 



~Gmi 



Ks 



where there are A'^vir particles in the halo, each of mass rrip. We 
sample TVfc particles from TVvir, \rij\ is the magnitude of the sep- 
aration between particles i at fl and j at fj, and the pref actor 
{N^ — TV)/(A| — Nk) accounts for particle sampling, e is the 
gravitational softening and Ks corresponds to the softening kernel 
used in GADGET2. For the Monte Carlo realisations in Figure[T]we 
set e to be vanishingly small, but for the simulations we use e as it 
is listed in table[T] 

Figure [T] shows \ W\ measured for Monte Carlo realisations 
of a halo with c =10 and rvir =200 kpc as a function of Avir- For 
comparison the horizontal dotted lines indicate the value of \W\ 
(±5%) we expect from equation l[6]l. If Avir ~ 300 or fewer, the 
measured \ W\ deviates from the expected \ W\ by greater than 5%; 
therefore we might regard TVcut = 300 as the lower limit on TVvir 
for a halo to be included in our sample. However, we adopt a more 
conservative TVcut = 600 in the remainder of this paper; this is be- 
cause the structure of simulated haloes ar e affected by finite grav- 
itational softening (cf. IPower et al.ll2003h . the y are are seldom (if 
ever) smooth and spherically symmetric (e.g. ISailin & Steinmet3 
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Figure 1. How many particles are required to measure accurately tlie 
gravitational binding energy of a dark matter halo?. Here we gener- 
ate Monte Carlo realisations of a NFW halo and calculate the gravitational 
potential energy of material within the virial radius. If there are too few 
particles within rvir, the potential energy will be inaccurate. 

l2005h , as we assumed in our simple cal culation, and there c an be a 
range of concentrations at a given mass teullock et al.l200lh . which 
will affect any estimate of W as inspection of equation ^ reveals. 
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Figure 2. Relationship between Virial Mass and Formation Redshift. 

Here we show how the formation redshift 2:torm varies with virial mass 
A/vir at z=Q. We determine Zfomi directly from a halo's merger tree - for a 
halo of A/vir identified at z=0, we identify the redshift Zfomj at which the 
mass of its main progenitor first exceeds half its virial mass at z=0. Data are 
binned using equally spaced bins in LogioA/vir^ filled circle and bars cor- 
respond to medians, upper and lower quartiles. The solid, upper and lower 
dashed curves corresponds to the median Zfomi and its up per and lower 
quarti les predicted by extended Press-Schechter theory (cf. iLacev & Cold 
1 19931) . 



3 QUANTIFYING MASS ASSEMBLY & MERGING 
HISTORY 

In this section we establish quantitative measures for a halo's mass 
accretion and merging histories, and we examine how these mea- 
sures relate to virial mass A/vir and formation redshift Zform- 

Quantifying Formation Redshift We begin our analysis by ver- 
ifying the correlation between virial mass A/vir and formation 
redshift Zform for ou r halo sample. We adopt the convention of 
ICole & LacevI j 19961) and define Zform as the redshift at which the 
mass of the main progenitor of a halo of mass Mvir{z) identified 
at z first exceeds Mvii[z)/2. This is equivalent to Zi/2. mb in the 
survey of halo formation redshift definitions examined bv lLi et al] 

feoosh . 

Our expectation is that more massive CDM haloes will assem- 
ble more of their mass at later times than their lower mass counter- 
parts and this is borne out by Figure[2] Here we show the variation 
of Zform with Mvir for our halo sample; the filled circles and bars 
indicate the medians and upper and lower quartiles respectively, 
within logarithmic mass bins of width 0.5 dex. The relationship be- 
tween the mean and median Zform with A/vir can be well approxi- 
mated by 



(Zfornr) ^ -0.22 log^o A/l2 + 1.06, 

and 



(8) 



Medzform ==; -0.23 logio Mi2 + 1.08, (9) 
where Mrz is Afvir in units of 10^^/i~^Mq. This is in very good 



agreement with the mean variation reported for the "Overall" sam- 
ple of halo es drawn from the M ille nnium and Millennium II sim u- 
lations (cf. lSpringel et al.ll2005l andlSoylan-Kolchin et al1l2009l re- 
spectively) in Table 3 of lMcBride et alj (l2009h . who found 



(■Sform) 



, Mi2 + 1.26. 



We show also the variation predicted by extended Press-Schechter 
(EPS) theory for our choice of CDM power spectrum - see the 
solid and dashed curves, indic ating the median, up per and lower 
quartiles of the distributions (cf. lLacev & ColJ 19931) . These curves 
were generated using realisations of 10^ Monte-Carlo merger trees 
for haloes with z=Q masses in the range 10^^° ^ MvIt / Mq ^ 
10^^'^. We note a slight but systematic offset between the medi- 
ans evaluated from the simulated haloes and those predicted by 
EPS theory, such that the simulated haloes tend to form earlier 
than predicted. This effec t has been reported prev iously by both 
Ivan den BoschI l l2002h and lMaulbetsch et alj l l2007h . 



Quantifying Recent Mass Accretion History Because more 
massive systems tend to form later than their less massive coun- 
terparts, it follows that the rate at which a halo assembles its mass 
should increase with increasing A^vir and decreasing Zform. The 
recent comprehensive study by iMcBride et al] ( |2009|) provides a 
useful fitting formula that captures the complexity of a halo's mass 
accretion history and allows haloes to be categorised into different 
Types I to IV, which depend on their growth rates. However, we 
adopt two simple well-defined measures of a halo's mass accretion 
rate that have a straightforward interpretation; 
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• (AAf/Af)At, the fraction of mass that has been accreted by a 
halo during a time interval At; and 

• a = 1/MdM/dz, the rate of fractional change in a halo's 
virial mass with respect to redshift over a redshift interval Az. 

Note that a i s equ ivalent to the a free parameter used in 
IWechsler et all ( |2002|) . We find that {AM/M)At and a are suf- 
ficient as simple measures of the mass accretion rate and we use 
them in the remainder of this paper. 

For the fiducial timescale At, we use twice the dynamical 
timescale Tdyn estimated at the virial radius, 



Tdyn 



-1/2 



70 



Gyrs 



(10) 



Note that Tdyn depends only on z and is the same for all haloes. 
For our adopted cosmological parameters, A{z) ~ 97 at z=0, 
and so At = 2 Tdyn — 5.6 Gyr which corresponds to a redshift 
interval of Az ~ 0.6 at z=0. Merging proceeds on a timescale 
Tmerge > Tdyn, with Tmcrgc — > Tdyn as the mass ratio of the merger 
decreases. Our adopted timescale of At = 2rdyn for the response 
of a halo to a merger is reasonable when compared to typical values 
of Tmeige / Td yn expected for haloes in cosm ological simulations, as 
estimated by'Sovla n-Kolchin et al] j2008lFI . 

We determine both (AM/M)Tdy„ and a directly from each 
halo's merger tree by tracking Mvir(«) of its main progenitor over 
the interval Az; a is obtained by taking the natural logarithm of the 
progenitor mass at each redshift and estimating its value by linear 
regression. Haloes that have high mass accretion rates will have 
{AM/M)At 1 and Of -oo. 

In Fig.[3]we show how a halo's mass accretion rate correlates 
with its virial mass and formation time. (AM /M)T^y^ (a) shows 
a steady monotonic increase (decrease) as Mvir increases over the 
range IO^^/i'^Mq < Myir < lO^'^h'^MQ. For example, inspec- 



0- 

reveals that < 15% of the virial mass of a 
10^^ Mq has been accreted since z ~ 0.5, 



tion of (AM/M)^^ 
halo with Mvi 

compared to ~ 50% for haloes with Mvir ~ 10^^/i"^Mq over 
the same period. (AA//M)Tjy„ (a) shows a similar increase (de- 
crease) with decreasing 2:form although it's interesting to note that 
the trend flattens off for haloes that form at 2 > 2. 

This analysis confirms our theoretical prejudice that more 
massive haloes and haloes that formed more recently tend to be 
the haloes with the measurably highest accretion rates. Reassur- 
ingly, our results are in good agreement wi th the findings of recent 
studies. For example. iMcBride et al.l | |2009|) examined the mass ac- 
cretion and merging histories of a much larger sample of haloes 
drawn from the Millennium and Millennium-II simulations and 
found that the mean instantaneous mass accretion rate varies with 
halo mass as compares fa vourably with our 

equivalent measure, (AM/M)^^^„ oc M°-3'^. lMaulbetschetal.l 
( 12007,) looked at halo accretion rates, normalised to their maximum 
masses, over the redshift interval 2=0.1 to for haloes with masses 
10^^ ^ Mvir/ft~^MQ ^ 10^'^ and found only a weak dependence 
on halo mass, with higher mass haloes have higher rates. This is 
consistent with with our results for q, whose median value changes 
by ~ 10% over the same range in halo mass. 



^ In particular, we refer to their equation 5 with values of j/jc{E) = 
0.5 and rc(E)/rvir that are consistent with the results of cosmological 
simulations. Here j is the specific angular momentum of a merging subhalo, 
jc{E) is the specific angular momentum of the circular orbit corresponding 
to the subhalo's orbital energy E, and tq^E) is the radius corresponding to 
this circular orbit. 
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Figure 3. Relationship between Recent Mass Accretion History, Halo 
Mass and Formation Redshift. For each halo of virial mass Mvir and 
formation redshift Zform identified at 2=0, we follow its merger tree back 
for one dynamical time Tdyn (— 4.6 Gyrs, A.z ~ 0.45) and characterise 
its mass accretion history using two measures. The first is (AAf/A/)Tjy„, 
the fraction of mass accreted over Tdyn (up per p anels), and the second is 
a, the average mass accretion rate of Wec hsler et al j i2002l) (lower panels). 
Data points and bars correspond to medians and upper and lower quartiles. 
Note that we use equally-spaced logarithinic bins in Mvir ™d Zform- 
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Figure 4. Frequency of Major Mergers and Dependence on Halo Mass 
and Formation Redshift. Here we determine the most significant merger 
of mass ratio <5max = Mace i^i ) / Mvir (2/ ) experienced by each halo since 
2=0.5, where Zi and Zf correspond to the initial and final redshifts. We then 
compute the fraction of haloes /(<5max) at a given virial mass (left hand 
panel) and given forination redshift (right hand panel) that have experienced 
mergers with mass ratios 5max in excess of 10% (filled circles), 20% (filled 
squares) and 50% (filled triangles). 
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Quantifying Recent Merger Activity Both (AAf /M )rdy„ and a 

provide useful insights into a halo's total mass accretion rate, but 
they cannot distinguish between smooth and clumpy accretion. In 
Fig. |4]we focus specifically on a halo's merger history by consid- 
ering the likelihood that a halo of a given A/vir (left hand panel) or 
zioTYn (right hand panel) has experienced at least one merger with 
a mass ratio 5max since 2=0.6. 

Each halo identified at z=Q has a unique merger history, which 
characterises not only how its Afvir grows as a function of time 
but also details of mergers it has experienced over time. Using 
this merger history, we construct the distribution of mass ratios 
of mergers 5 experienced by a halo of a given Afvir or Zform be- 
tween ^ 2: < 0.6. We define 5 = Mi^^ccizi, z/)/Afvir(z/) where 
Afacc(2i, z/) is the mass of the less massive halo prior to its merg- 
ing with the more massive halo, Afvir(z/) is the virial mass of the 
more massive halo once the less massive halo has merged with it, 
and Zi > Zf and Zf are the redshifts of consecutive simulation 
snapshots. The maximum value of 5 for a given halo gives us its 
^max and we use this to compute the fraction of haloes of a given 
Afvir or Zform that have 5max in excess of 10% (filled circles), 20% 
(filled squares) and 50% (filled triangles). 

Fig. m reveals that mergers with higher mass ratios (i.e. mi- 
nor mergers) are more common than mergers with lower mass ra- 
tios (i.e. major mergers), independent of M^ir and Zform, and that 
more massive (older) haloes tend to experience more mergers than 
their lower mass (younger) counterparts. For example, the likeli- 
hood that a IO^^/i'^Mq galaxy-mass halo experiences a merger 
with 5niax > 10% is ~ 35%, compared to 25%(10%) for (5max > 
20%(50%). In contrast, the likelihood that a 10"/i"^Mq cluster- 
mass halo experiences mergers with 5max > 10%(20%, 50%) 
is ~ 60% (40%, 20%). Interestingly, we find that the fraction of 
haloes that have experienced a merger more significant than 5max 
increases with Afvir approximately as /((5max) oc M^i^^ . 

These res u lts are broadly in agreement with the findings of 
iFakhouri et al.l ( l2010h . Inspection of the leftmost panel of their 
Figure 7 shows the mean number of mergers with mass ratios 
greater than 1:10 and 1:3 between z=0 and z ~ 0.6 increases 
with increasing halo mass, such that a 10^^ (10^*)/i~'^Mq has 
a likelihood of ~ 40% (~ 80%) to have experienced a merger 
with 5niax > 10%, and a likelihood of ~ 20% (~ 40%) to have 
experienced a merger with (5max > 33%. 

In Fig. |5]we show the full (cumulative) distributions of (5max 
for haloes split into bins according to Afvir (left hand panel) 
and 2:form (right hand panel); note that we consider only haloes 
with (5max 10%. Interestingly this figure reveals that the 
probability distribution of 5max is insensitive to Afvir, but depends 
strongly on jzform. For example, the median (Jmax.mod — 0.3, 
independent of Afvir whereas it increases from 5max,med — 0.2 
for haloes with Zform 0.5 to 5max,med — 0.4 for haloes 
with 0.25 ^ zform ^ 0.5 and (Jmax.mod — 0.7 for haloes with 

Zform 0.25. 

Figs. [3] to |5] demonstrate that there is a strong correlation at 
z=0 between a halo's virial mass Afvir, its formation redshift Zform 
and the rate at which it has assembled its mass through accretion 
and merging over the last rdyn or equivalently Az ~ 0.6. We 
use these results in !j5] where we investigate the degree to which 
a halo's Afvir, Zform and mass accretion rate affect the degree to 
which it is in dynamical equilibrium. 
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Figure 5. Cumulative Distribution of 5max as a Function of Halo Mass 
and Formation Redshift. We show how the fraction of haloes whose most 
significant merger's mass ratio is less than (5,nax, as a function of virial mass 
(left hand panel) and formation redshift (right hand panel). Note that we 
select only haloes that have <5niax ^ 10%, and we consider only mergers 
between 2=0.5 and 2=0. In the key, the numbers in brackets correspond to 
the lower and upper bounds in A/vir and 2form ■ 



4 QUANTIFYING DYNAMICAL EQUILIBRIUM 

In this section we describe the two commonly used quantitative 
measures for a halo's dynamical state, the virial ratio rj and the 
centre of mass offset Ar, and we examine their relationship with 
virial mass Afvir and formation redshift Zform. 



4.1 The Virial Ratio 77 

The virial ratio rj is commonly used in cosmological A^-body simu- 
lation s as a measure of a halo's dynamical state (e.g.lCole & Lace 



1 9961: iBett et al 



I200I lNetoet al.1 12007| : iKnebe & Power! 12001 



Davis et al.ll201 Ih . It derives from the virial theorem. 



i^ = 2r + iy + i5„ 



(11) 



where I is the moment of inertia, T is the kinetic energy, W = 
Ei^.r is the virial, and Es is the surface pressure integrated over 
the bounding surface of the volum e within which I, T and W are 
evaluated (cf. lChandrasekhail I96lh . Provided the system is isolated 
and bounded, the virial W is equivalent to the gravitational poten- 
tial energy. While not strictly true for haloes that form in cosmo- 
logical A'^-body simulations, the convention has been to evaluate W 
as the gravitational potential energy with this caveat in mind (e.g. 
ICole & Lacevlll996h . We follow this convention and treat W as the 
gravitational potential energy computed using equation 0. 

If the system is in a steady state and in the absence of sur- 
face pressure, equation Jilt reduces to 2 T + = 0, which can 
be wr itten more compactly as 2T/\W\ = 1 (e.g. ICole & LacevI 
We refer to the ratio rj — 2 T/\W\ as the virial ratio and we 
expect ry — >■ 1 for dynamically relaxed haloes. However, we might 
expect Es to be important for haloes that form in cosmological 
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Figure 6. Correlation between Virial Ratios r? and r;' . We bin all haloes 
in our sample at z=0 according to their r] and evaluate the median rj' within 
each bin. The upper and lower quaililes of the distributions in r) and r]' are 
indicated by bars. 



Figure 7. Correlation between Centre-of-Mass Offset Ar and Virial 
Ratios ri and rj' . We can clearly see the relation which is confirmed by 
measuring a Spearman rank coefficient of 0.45 whereas we find an anti- 
correlation with Spearman rank coefficient of -0. 1 8 for rj' . 



4.2 The Centre-of-Mass Offset Ar 



A'^-body simulations; in this case lShaw et aD (l2006h have proposed 
modifying the virial ratio to obtain 



r^' = {2T-E,)/\W\ 



(12) 



We calcul ate both T and W using all material within r^ir, while 
we follow IShaw et alj j2006l) by computing the surface pressure 
contribution from all particles that lie in a spherical shell with inner 
and outer radii of 0.8 and 1.0 rvi,-, 



Ps = —T.^{miVi) 



(13) 



here V corresponds to the volume of this shell and Vi are the 
particle velocities relative to the centre of mass velocity of the 
halo. The energy associated with the surface pressure is therefore 

Es ■ 

Figure[6]shows how the median 77 and 77' for the haloes in our 
sample compare, with bars indicating the upper and lower quartiles 
of the distributions. We might expect that rj' ~ 1 and insensitive 
to variation in -q; however, this figure reveals that the relationship 
between rj and rj' is not so straightforward. Haloes that we would 
expect to be dynamically relaxed, with 7; ~ 1, have values of rj'iO, 
suggesting that Es t ends to over-correct . Si milar behaviour has 
been noted in both Knebe & Poweii l2008h and lDavis etal] i201lh 
for high redshift haloes (z > 1). The relation between the median 
7] and rj' is flat rj < 1.25 but rises sharply from rj' ~ 0.9 to peak 
at rj' ~ 1.05 before declining sharply for rj > 1.4 to a median of 
rj' ~ 0.8 in the last plotted bin. Interestingly, the width of the rj' 
distribution increases with rj; if rj tracks recent major merging ac- 
tivity as we expect, then this suggests that rj' - and consequently 
the surface pressure correction term Es - is sensitive to mergers 
but in a non-trivial way. 



Another commonly used measure of a halo's dynamical state is the 
centre-of-mass offset Ar-, 



Ar = 



(14) 



which measures the separation between a halo's centre-of-density 
rccn (calculated as described in ^2.2\ and its centre-of-mass, 
calculated using a l l material within ryj r , normalised by ryir (cf . 
Crone et al.l 1 19961: iThomas et alj 1 19981 | 200ll: iNeto etai] |2007|; 



Maccio et al 



20071 : iD'Onghia & Navarroll2007l) . Ar is used as a 



substructure statistic, providing an estimate of a halo's deviations 
from smoothness and spherical symmetry. The expectation is 
that the smaller t he Ar, the more relaxed the halo; for example, 
iNeto et al.l ( l2007h define dy namically relaxed halo e s to b e those 
with Ar ^ 0.07, while iD'Onghia & Navarrcl ^OOf) adopt 
Ar ^ 0.1. Maccio et alj ( l2007h favoured a more conservative 
Ar ^ 0.04 based on a thorough analysis. 

We can get a sense of how well Ar measures the dynamical 
state of a halo by comparing it to rj and rj' . In Figure |7] we plot the 
median rj and rj' (filled circles and squares respectively) against the 
median Ar ; as before, bars indicate the upper and lower quartiles of 
the distributions. This figure shows that both rj and rj' correlate with 
Ar - but in different senses; as Ar increases, rj increases while rj' 
decreases. The increase (decrease) is a gradual one; for example, 
for Ar < 0.04, the median rj is flat with a value of ~ 1.05, but for 
Ar > 0.04 there is a sharp increase and Ar > 0.1, rj ~ 1.2. Al- 
though direct comparison is difficult, a similar trend can be gleaned 
from Figure 2 of Neto et al. (200?,). We use the Spearman rank cor- 
relation coefficien t to assess the strength of t he correlation between 
Ar and 77(77') (cf. iKendall & Gibbonsll 1990 1. and find strong pos- 
itive and negative correlations for 77 (Spearman rank coeffecient 
r=0.97) and 77' (r=-0.95) respectively. 
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Figure 8. Relationship between Dynamical State and Halo Mass and 
Formation Redshift. For haloes identified at 2=0, we plot the median rj and 
r)' versus M^ir (left hand panel) and Zform (right hand panel) using equally 
spaced bins in LogioMvir and Zform- Data points and bars correspond to 
medians and upper and lower quartiles. 
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Figure 9. Relationship between Centre-of-Mass Offset, Halo Mass and 
Formation Redshift. For haloes identified at 2=0, we plot the median 
centre-of-mass offset Ar versus Mvir (left hand panel) and Zf^iui (right 
hand panel) using equally spaced bins in LogioMvir and Zfoj-ni- Data 
points and bars correspond to medians and upper and lower quartiles. 



This is suggestive - as we show below, Ar correlates more 
strongly with merging activity than either of r; or r;' (cf. Figure [T3l 
in Sj5](. Both Ar and iq increase with strength of merging activity, 
whereas 77' appears to be over-corrected by Es (as we have noted 
above). From this we conclude that Es (as we evaluate it) correlates 
with significant merger activity, which is confirmed by a Spearman 
rank coefficient o f 0.38 for the corre lation between Es and (5max- 

Interestingly lOavis et alj j201 ih examined the correlation be- 
tween Ar and 77' for high redshift haloes (z > 6) and noted a 
tendency for haloes with small values of ry' to have larger values of 
Ar. Inspection of their Figure 4 shows that this is true for haloes 
with 0.4 < Ar< 10; for Ar < 0.4 the relation with 77' is flat. 
lDavisetal.U201lh argue that, because there is no systematic shift 
in r]' for Ar < 0.1, Ar is not a useful measure of dynamical state 
at high redshifts. 



4.3 Dependence of Dynamical State on Afvii and zform 

In Figures [8] and [9] we examine how r],rj' and Ar vary with i\fvir 
(left hand panels) and Zform (right hand panels) for the halo pop- 
ulation at z=Q. Haloes are sorted in bins of equal width in mass 
(A logio M=0.5 dex) and redshift (Az=0.25), and we plot the me- 
dian r//r/'/Ar within each bin against the median A/vir/zform; bars 
indicate the upper and lower quartiles of the respective distribu- 
tions. For reference, we also plot a horizontal dotted line in each 
panel of Figure[8]to indicate a virial ratio of unity. 

Because more massive haloes tend to form at later times, and 
because these haloes tend to assemble a larger fraction of their mass 
more recently, we expect that 77 and Ar should increase with in- 
creasing Mvir and decreasing Zform- This is borne out in Figures [8] 
an d|9] We find that the mean and median Ar increases steadily with 
increasing Mvir as 
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(logio Ar) = -1.47 + O.OSlogio M12 



(15) 



Figure 10. Distribution of virial ratios. Here we show the correlation be- 
tween halo mass and virial ratio rj and rj' at redshift 2 = 0. 



and 

Med login Ar = -1-49 0.09 logjo M12 (16) 

where, as before, M12 is Mvir in units of 10^^ h~^Mf^. This is 
consistent with the result of iThomas et alj ( 1200 iL see their Figure 
9), who found a similar trend for Ar to increase with Migo for a 
sample of cluster mass haloes (lO" < Miso/{h~^MQ) < 10^^) 
in a rCDM model. Their typical values of Ar are offset to higher 
values than we find, but this can be understood as an effect of A, 
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the merging rate being suppressed in the ACDM model compared 
to the rCDM model. Similarly, Ar varies strongly with Zfam', for 
Zform > 1 we find that Ar oc (1 + compared to Ar oc 

(l + z)-0-^forzfo,n, > 1. 

The mean and median rj exhibit similar behaviour, increasing 
with increasing Afvir, albeit weakly, as. 



(logio r?) = 0.05 + 0.016 logio M12 
and 

Med logio -n = 0.04 + 0.019 logio Mi 



(17) 



(18) 



This means that rj is systematically greater than unity for all Mvir 
that we consider - ~ 1.15 for a typical 10^^/i~^Mq halo, com- 
pared to 77 ~ 1.25 for a typical 10^^/i~^Mq halo. The same grad- 
ual increase in with decreasing Zform is also apparent. 

As we might have anticipated from inspection of Figures[6]and 
|7] 77' is systematically smaller than unity. Its variation with Mvir 
is negligible (oc M^i^'^"'^; a littl e surprising, when com pared to 
oc Mvif^^ at 2:=1, as reported bv lKnebe & Powej[20o3) but there 
is a trend for the median ri' to decrease with decreasing Zform. This 
makes sense because Es increases with the significance of recent 
mergers and haloes that have had recent major mergers tend to have 
smaller Zfovm- This effect is also noticeable in the width of the ri' 
distributions in each bin (as measured by the bars), which are larger 
than than the corresponding widths of the rj distribution. 

We look at this effect in more detail by plotting the distribu- 
tions of rj and rj' shown in Figure [TO] Here it is readily apparent 
that there is a systematic shift towards larger rj as A/vir increases. 
Interestingly the rj' distribution remains centred on r;' ~ 0.9, but 
it spreads with increasing A/vir; again, this suggests the sensitivity 
of ri' to recent merging activity. 

Figs. |6] to [To] demonstrate that there is a strong correlation 
at z=0 between a halo's virial mass A/vir, its formation redshift 
-Zform and its dynamical state, as measured by the virial ratio ri and 
the centre-of-mass offset Ar. In contrast, the correlation with ri' 
is more difficult to interpret, especially when 77 is large. In these 
cases, we expect significant merging activity and as we note above, 
the correction by the surface pressure term Es increases the width 
of the original i] distribution by a factor of ~ 2-3. It's also note- 
worthy that the median rj' is systematically offset below unity. For 
this reason we argue that 77' is not as useful a measure of a halo's 
dynamical state as rj. 



5 LINKING MASS ASSEMBLY AND DYNAMICAL 
STATE 

We have established quantitative measures of a halo's mass assem- 
bly and merging history and its dynamical state in the previous two 
sections, and we have investigated how these relate separately to 
a halo's virial mass Afvir and its formation redshift Zform. In this 
final section we examine the relationship between a halo's mass 
assembly history and its dynamical state directly. 

In Figures [TT] and [T2] we show explicitly how a halo's re- 
cent mass accretion and merging history impacts on its virial ra- 
tio. As in section[3l we quantify a halo's mass accretion history by 
{AM/M)r^ „, the fractional increase in a halo's mass over the pe- 
riod Tdyn (equivalent to a redshift interval Az ~ 0.6 at 2=0), and 
Q, the mean accretion rate over the period Tdyn). We use 5max, the 
mass ratio of the most significant merger experienced by the halo 
over Tdyn, to characterise a halo's recent merging history. 
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Figure 11. Relationship between Recent Mass Accretion and the Virial 
Ratio. Here we investigate how (AM/M)rdy„. the fraction of mass ac- 
creted over Tdyn (left hand panel), and a, the mean accretion rate over 
Tdyn (right hand panel), correlate with the standard (r], filled circles) and 
corrected (t)', filled squares) virial ratio. Data points correspond to medians 
and bars correspond to the upper and lower quartiles. 
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Figure 12. Relationship between Most Significant Recent Merger and 
the Virial Ratio. Here we investigate how (5max, which measures the mass 
ratio of the most significant recent merger since z=0.5, correlates with 
the standard (rj) and corrected (rj') virial ratios respectively. Filled circles 
(squares) con'espond to medians oirj (r)' ), while bars indicate the upper and 
lower quaililes. 
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Figure 13. Correlation between centre-of-mass offset Ar and Recent 
Merging and Accretion History. Here we examine whether Ar correlates 
with the fraction of inass accreted over T^y^, {AM/M)t^^^, the mean 
accretion rate a and the most significant inerger (Smax- Filled circles corre- 
spond to medians in the respective bins; bars represent the upper and lower 
quartiles of the distribution. 



Figure 14. Relationship between rj, Ar and zg^^^, the redshift of the 
most significant recent major merger. We identify all haloes in our sam- 
ple at z=0 with 5niax JSJ 1/3 and identify the redshift zs^_^^^ at which (5max 
occurred. Both rj and Ar are evaluated at z=0. Filled circles and bars cor- 
respond to medians and upper and lower quartiles. 



We expect that the standard virial ratio rj should increase with 
increasing mass accretion rate and decreasing mass ratio of most 
significant merger, which is in good agreement with the behaviour 
that we observe. In particular, the median variation of rj with a 
and (5max can be well approximated by logj^Q ~ 0.004 — 0.126a 
and rj ~ 1.251^1-^; the corresponding variation of Ar can be well 
approximated by 0.01 — 0.1a and 0.1(5,nax- 

Interestingly we note that the median corrected virial ratio ri' 
declines with increasing mass accretion rate and mass ratio of most 
significant merger. Both correlations indicate that merger events 
lead to a state that is less virialised, but, as we have noted al- 
ready, the inclusion of the surface pressure term over-corrects the 
virial ratio. We see in Figure [m that for (AM/M).r<iy„ < 0.2, 
both the median rj and rj' are flat; rj ~ 1.05 whereas the median 
rj' ~ 0.85. Above {AM/M)-r^^^ ~ 0.2, the median increases 
sharply whereas it is the width of the r]' distribution that shows the 
sharp increase. Comparison with Figure[72lprovides further insight 
- the median rj (rj') shows a gradual increase (decrease) with in- 
creasing 5max, starting at 77 ~ 1.05 (?/ ~ 0.9) for 5max ~ 0.02. 
For (5max=0.1, rj ~ 1.1 (rj' ~ 0.85). However, whereas the width 
of the rj distribution is largely insensitive to 5max, the width of the 
rj' distribution increases rapidly, bearing out our observations in the 
previous section. 

In Figure[T3]we show how Ar varies with {AM/M)t^ „, a 
and (5niax. This reveals that Ar increases with increasing mass ac- 
cretion rate and mass ratio of most significant recent merger, as we 
would expect. Although the scatter in the distribution is large, we 
can identify the remnants of recent major mergers (Jmax ^ 30%) 
as haloes with Ar > 0.06. Haloes that have had relatively qui- 
escent recent mass accretion histories ((AM/M)Tdy„ J; 0.1, 
5max < 10%) have Ar < 0.04. 



Merging Timescale & Dynamical State We conclude our analy- 
sis by investigating the timescale over which the effect of a merger 
can be observed in the virial ratio 77 and the centre-of-mass offset 
Ar. 

In Figure [14] we investigate how a typical halo's 77 (upper 
panel) and Ar (lower panel), measured at 2 = 0, correlates with 
the redshift at which the halo suffered it's most significant merger, 
-^■Smax ■ For clarity, we focus on haloes for which (5max > 1/3, al- 
though we have verified that our results are not sensitive to the 
precise value of (5max that we adopt; filled circles correspond to 
medians and bars indicate upper and lower quartiles. The median 
rj increases with decreasing zs^^^ for zg^^^ > 1 before peak- 
ing at zg^^^ ~ 0.3 and declining at lower zg^^^. The median 
Ar shows a similar steady increase with decreasing zg^^^^ below 
•^4max ~ 0-8 although there is evidence that it peaks at zg^.^^ ~ 0.4 
before declining at lower zg^^^. The redshift interval correspond- 
ing to zg^^^^ ~ 0.3 represents a time interval of At ~ 4.3 Gyrs or 
~ 1.5rdyn. 

This is consistent with the finding of iTormen et all ( Il997h . 
who examined the velocity dispersion Vims of material within rvir 
of simulated galaxy cluster haloes (see their Figure 5). They noted 
that merging leads to an increase in Wrms of the main (host) halo 
because the merging sub-halo acquires kinetic energy as it falls in 
the potential well of the more massive main halo. The peak in Wrms 
corresponds to the first pericentric passage of the subhalo, after 
which Urms dcclines because subsequent passages are damped, 
and so the main halo relaxes. This will occur on a timescale of 
order ~ 1 — 2rdyn, which is consistent with the peak in rj at 
zg^.^^ ~ 0.3. We would expect to see a peak in Ar on roughly the 
merging timescal e Tmcrgc, which as we noted in §[3]is comparable 
to 1 - 2 Tdyn (cf. lBovlan-Kolchin et al.ll2008l) . 

We can take this analysis a little further by looking at the de- 
tailed evolution of rj and Ar over time. In Figure[T5]we plot the red- 
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Figure 15. Response of M^ir, V and Ar to a major merger. We include 
all haloes with <5max 1/3 at z=0 and plot the redshift variation of Mvir, 
7? and Ar against the time since the major merger, normalised by the dy- 
namical time of the halo at the redshift at which the merger occuiTed. Filled 
circles and bars correspond to medians and upper and lower quartiles. while 
curves correspond to the histories of 5 individual haloes. 



shift variation of Mvir (normalised to its value at z=0; lower panel), 
rj (middle panel) and Ar (upper panel) against the time since ma- 
jor merger, normalised by the dynamical time Tdyn estimated at the 
redshift at which the merger occurred, zg^^^ . Medians and upper 
and lower quartiles are indicated by filled circles and bars. For il- 
lustrative purposes, we show also the redshift variation of Afvir, V 
and Ar for a small subset of our halo sample (red, blue, green, cyan 
and magenta curves). As in Figure[T4l we adopt (5max Si 1/3. 

Our naive expectation is that both rj and Ar should increase 
in response to the merger, peak after At ~ Tdyn and then re- 
turn to their pre-merger values. If this behaviour is typical, then 
we expect pronounced peaks in the median values of rj and Ar at 
A t/rdyn — 1- However, it is evident from Figure [T5] that there is 
no significant difference between the median rj and Ar pre- and 
post-major merger, and so our naive expectation is not borne out by 
our results. 

This is not surprising if on e inspects histo ries for rj and Ar for 
individual haloes, in the spirit of lTormen et al.l (IT997): i] and Ar in- 
crease following a major merger, but the behaviour is noisy (reflect- 
ing e.g. differences in orbital parameters of merging subhaloes, the 
redshift dependent virial radius, dependence on environment, etc..) 
and the timescale of the response varies from halo to halo - simply 
averaging or taking the median washes any signal away. Neverthe- 
less it is worth looking at this in more detail, which we shall do in 
a forthcoming paper 



6 SUMMARY 

The aim of this paper has been to quantify the impact of a dark mat- 
ter halo's mass accretion and merging history on two measures of 
dynamical state that are commonly used in cosmol ogical A^-body 
simul ations, namely the virial ratio 77 = 2T/\W\ (cf. lCole & LacevI 

Il996h and the centre-of-mass offset — 1 rccn — fcm l/rvir (cf. 



ICrone et al.ll996l : lThomas et aljl998l,l200lh . The virial ratio 77 de- 
rives from the virial theorem and and the expectation is that rj ^ 1 
for dynamically equilibrated haloes. The centre- of-mass offse t Ar 
can be regarded as a substructure statistic ( Tho mas et alj20oi1) that 
provides a convenient measure of how strongly a halo deviates from 
smoothness and spherical symmetry. We expect a halo's dynam- 
ical state and its mass assembly history to correlate closely. Un- 
derstanding how is important because the degree to which a halo 
is dynamically equilibrated affects the reliability with which we 
can measure characteristic propert ies of its structure, such as the 
concentration parameter Cvtr (e.g. iMaccio et al.ll2007l : iNeto et al.l 
[2007;Prada etal.20li^ and kinematics, such as the spin parameter 
A (e.g . IGardneril200ll : lD'Onghia & Navarrdl2007l ; iKnebe & Powe3 
I2OO8I) . For this reason, it is desirable to establish quantitatively how 
well they correlate and to assess how Tj and Ar can help us to char- 
acterise how quiescent or violent a halo's recent mass assembly 
history has been. 

Our key results are that rj and Ar show strong positive corre- 
lations with each other (cf. Figure|7]( - as 77 increases for a halo, so 
too does Ar ~ and that both are useful indicators of a halo's mass 
recent mass accretion and merging history. For example, 77 and Ar 
correlate strongly with (5max, which measures the significance of a 
halo's recent merging activity; haloes with rj < 1.1 (cf. Figure [T2]l 
and Ar < 0.04 (cf. Figure [T3t will have quiescent recent mass as- 
sembly histories - they are unlikely to have suffered mergers with 
mass ratios greater than 1:10 over the last few dynamical times. 

In contrast, interpreting the corrected virial ratio rj'={2T — 
Es)/\W\, where Es is the surface pressure energy, is problem- 
atic (at least insofar as we have impleme nted it here, which fol- 
lows the pr escription of IShaw et alj|20 06' and has been applied in 
iKnebe &"P ower 2008 and lOavis et al. [201 1). In principle, rj' should 
account for the approximation that is made when we define a halo 
to be a spherical overdensity of Avir times the critical density at 
a particular redshift. As we noted in § |2] haloes are more complex 
structures than this simple working definition gives them credit for, 
and by defining the halo's extent by the virial radius r^ii the likeli- 
hood is that material that belongs to the halo will be neglected. By 
correcting the virial ratio 77 for what is effectively a truncation of the 
true halo, the corrected virial ratio rj' takes account of the "missing" 
kinetic energy. However, our results imply that the correction itself 
(the surface pressure energy Es) is sensitive to a halo's merging 
history, and that it increases with increasing (5max (cf. Figure [T2J. 
For this reason we would caution against the use of 77' to identify 
dynamically relaxed haloes, at least in the form that is currently 
used. 

Interestingly, we find that systems with violent recent mass 
assembly histories (most significant merger with a mass ratio 
5max ^ 1/3 between < 2 < 1) have values of 77 and Ar 
(as measured at z=0) that peak at zs^^^ — 0.3 — 0.4, which 
corresponds to a timescale of ~ 1.5 Tdyn (cf. Figure II 4t. Thi s 
is consistent with the earlier analysis of iTormen et alj ( 1 19971) . 
who found that the velocity dispersion v^rna of material within 
the virial radius - which is linked to the virial ratio 77 ~ peaks 
on first closest approach of the merging sub-halo with the centre 
of the more massive host halo. This should occur on a timescale 
of ~ 1 — 2 Tdyn, after which Umis and 77 should dampen away. 
Similar arguments can be made for Ar. We note that these 
arguments can be made in a statistical sense, but if we look 
at the merging histories of individual haloes, the behaviour of 
?7 and Ar is much more complex, and as we demonstrate a 
simple timescale for their response to a major merger is difficult 
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to define (cf. FigurefTSll. We shall return to this topic in future work. 

What is the significance of these results? Structure formation 
proceeds hierarchically in the CDM model and so we expect to 
find correlations between virial mass Mvir and formation redshift 
2form (cf. Figure[2}, which in turn result in positive correlations be- 
tween Afvir/2:form and ri/Ar. (cf. Figures[8]and|9}. This means that 
more massive haloes and those that formed more recently are also 
those that are least dynamically equilibrated, a fact that we should 
be mindful of when characterising the halo mass dependence of 
halo properties that are sensitive to dynamical state (e.g. Cvir and 
A). It's worth noting that the correlation between Mvir and rj is 
stronger than the correlation between Alvh and Ar; the median 77 
rises sharply with A/vir and there is no overlap between the width 
of the distributions of r; is the lowest and highest mass bins. In con- 
trast, the median Ar in the highest mass bin lies in the high- Ar tail 
of the lowest mass bin. 

This is interesting because rj as it is usually calculated 
depends on W, which is sensitive to the precise boundary of the 
halo. Correcting for the surface pressure term does not appear 
to help, as we point out ~ indeed, the surface pressure term 
itself correlates with merging activity. This points towards an 
ambiguity in the use of - as we note, it rarely if ever satisfies 
ri=l. We discuss this point in a forthcoming paper, but we note 
that even i n ideal situations, what one computes for ? ? depends 
on rvir (cf. lCole & Lacevlll996lJLokas & Mamoij|200ll) - and so 
applying a flat cut based on a threshold in 77 alone risks omitting 
massive haloes that might otherwise be considered dynamically 
equilibrated. For this reason we advocate the use of Ar in cosmo- 
logical A'^-body simulations as a more robust measure of a halo's 
dynamical state; its calculation is computationally inexpensive, it 
is well defined as a quantity to measure, and its interpretation is 
both clear and straightforward. We find that Ar < 0.04, which 
corresponds to a (5max 0.1, should be sufficient to pick out 
the most dynamically relaxed haloes in a simulation volume at z=0. 

Although our focus has been fixed firmly on haloes in cos- 
mological simulations, we note that our results have observational 
implications. Whether or not an observed system - for example, a 
galaxy cluster - is in dynamical equilibrium will affect estimates 
of its dynamical mass if we assume a luminous trace r population 
that is in dynamical equilibrium (e.g. lPiffaretti & 'V^darnini 20()8h . 
Similarly, studies that seek to reconstruct a galaxy cluster's re- 
cent merg ing history tend to use signatures of disequilibrium (e.g. 
ICassano e t al. 2010). The most obvious measure of disequilibrium 
is the centre of mass offset Ar, or its projected variant. Although 
a more careful study in which we mock observe our haloes (and a 
seeded galaxy population) is needed, our results suggest that Ar 
could be used to infer the redshift of the last major merger (cf. Fig- 
ures[T3]and Figure|9l although care must be taken as Figure [Tslre- 
veals). Observationally, this would require measurement of, for ex- 
ample, projected displacements between gas and dark matter from 
gravit ational tensing and X-ray studies. We note that IPoole et al.l 
( 12006 ) have already tested this idea using idealised hydrodynamical 
simulations of mergers between galaxy clusters and found that the 
centroid offset between X-ray and projected mass maps captures 
the dynamical state of galaxy clusters well, but it is interesting to 
extend this idea using cosmological hydrodynamical simulations of 
galaxy groups and clusters. This will form the basis of future work. 



ACKNOWLEDGMENTS 

CP acknowledges the support of the STFC theoretical astrophysics 
rolling grant at the University of Leicester. AK is supported by the 
Spanish Ministerio de Ciencia e Innovacion (MICINN) in Spain 
through the Ramon y Cajal programme as well as the grants AYA 
2009-13875-C03-02, AYA2009-12792-C03-03, CSD2009-00064, 
and CAM S2009/ESP-1496. He further thanks the Aluminum 
Group for chocolates. SRK acknowledges financial support from 
Swinburne University of Technology's Centre for Astrophysics and 
Supercomputing's visitor programme. He acknowledges support 
by the MICINN under the Consolider-Ingenio, SyeC project CSD- 
2007-00050. The simulations presented in this paper were carried 
out on the Swinburne Supercomputer at the Centre for Astrophysics 
& Supercomputing, the Sanssouci cluster at the AIP and the ALICE 
supercomputer at the University of Leicester. 



This paper has been typeset from a TgX/ UrgX file prepared by the 
author. 



REFERENCES 

Bailin J., Steinmetz M., 2005, Apl, 627, 647 

Bett P, Eke V., Frenk C. S., Jenkins A., Helly J., Navarro J., 2007, 

MNRAS, 376, 215 
Boylan-Kolchin M., Ma C.-P, Quataert E., 2008, MNRAS, 383, 

93 

Boylan-Kolchin M., Springel V., White S. D. M., Jenkins A., 

Lemson G., 2009, MNRAS, 398, 1150 
Bullock J. S., Kolatt T. S., Sigad Y., Somerville R. S., Kravtsov 

A. v., Klypin A. A., Primack J. R., Dekel A., 2001, MNRAS, 

321,559 

Cassano R., Ettori S., Giacintucci S., Brunetti G., Markevitch M., 

Venturi T., Gitti M., 2010, ApJ, 721, L82 
Chandrasekhar S., 1961, Hydrodynamic and hydromagnetic sta- 
bility. International Series of Monographs on Physics, Oxford: 
Clarendon, 1961 
Cole S., Lacey C, 1996, MNRAS, 281, 716 
Crone M. M., Evrard A. E., Richstone D. O., 1996, ApJ, 467, 489 
Davis A. J., D'Aloisio A., Natarajan P, 201 1, MNRAS, pp 950-+ 
Diemand J., Kuhlen M., Madau P, 2007, ApJ, 667, 859 
Diemand J., Kuhlen M., Madau P., Zemp M., Moore B., Potter D., 

Stadel J., 2008, Nature, 454, 735 
Diemand J., Moore B., 2009, Arxiv e-prints ( |arXiv:0906.'4340t 
D'Onghia E., Navarro J. F, 2007, MNRAS, 380, L58 
Eke V. R., Cole S., Frenk C. S., Patrick Henry J., 1998, MNRAS, 
298, 1145 

Fakhouri O., Ma C.-P, 2008, MNRAS, 386, 577 
Fakhouri O., Ma C, 2010, MNRAS, 401, 2245 
Fakhouri O., Ma C.-P, Boylan-Kolchin M., 2010, MNRAS, 406, 
2267 

Gao L., Navarro J. F, Cole S., Frenk C, White S. D. M., Springel 
v., Jenkins A., Neto A. F, 2008, MNRAS, 387, 536 

Gao L., Frenk C. S., Boylan-Kolchin M., Jenkins A., Springel V., 
White S. D. M., 2011, MNRAS, 410, 2309 

Gardner J. P, 2001, ApJ, 557, 616 

Gin S. P D., Knebe A., Gibson B. K., 2004, MNRAS, 351, 399 
Hetznecker H., Burkert A., 2006, MNRAS, 370, 1905 



Kendall M., Gibbons J. D., 1990, Rank Correlation Methods, 5 

edn. A Charles Griffin Title 
Knebe A., Power C, 2008, ApJ, 678, 621 
Knollmann S. R., Knebe A., 2009, ApJS, 182, 608 
Lacey C, Cole S., 1993, MNRAS, 262, 627 
Li Y., Mo H. J., Gao L., 2008, MNRAS, 389, 1419 
Lokas E. L., Mamon G. A., 2001, MNRAS, 321, 155 
Maccio A. V., Dutton A. A., van den Bosch F. C., Moore B., Potter 

D., Stadel J., 2007, MNRAS, 378, 55 
Mandelbaum R., Seljak U., Hirata C. M., 2008, Journal of Cos- 
mology and Astroparticle Physics, 8, 6 
Maulbetsch C, Avila-Reese V., CoHn P., Gottlober S., Khalatyan 

A., Steinmetz M., 2007, ApJ, 654, 53 
McBride J., Fakhouri O., Ma C.-R, 2009, MNRAS, 398, 1858 
Navarro J. F, Frenk C. S., White S. D. M., 1997, ApJ, 490, 493 
Navarro J. F, Ludlow A., Springel V., Wang J., Vogelsberger M., 
White S. D. M., Jenkins A., Frenk C. S., Helmi A., 2010, MN- 
RAS, 402, 21 

Neto A. F, Gao L., Bett R, Cole S., Navarro J. F, Frenk C. S., 
White S. D. M., Springel V., Jenkins A., 2007, MNRAS, pp 922- 
-I- 

Piffaretti R., Valdamini R., 2008, A&A, 491, 71 

Poole G. B., Fardal M. A., Babul A., McCarthy I. G., Quinn T., 

Wadsley J., 2006, MNRAS, 373, 881 
Power C, Navarro J. F, Jenkins A., Frenk C. S., White S. D. M., 

Springel V., Stadel J., Quinn T., 2003, MNRAS, 338, 14 
Prada F, Klypin A. A., Cuesta A. J., Betancort-Rijo J. E., Primack 

J., 201 1. ArXiv e-prints 
Romano-Diaz E., Hoffman Y., Heller C, Faltenbacher A., Jones 

D., Shlosman L, 2007, ApJ, 657, 56 
Seljak U., Zaldarriaga M., 1996, ApJ, 469, 437 
Shaw L. D., Weller J., Ostriker J. P, Bode R, 2006, ApJ, 646, 815 
Springel V., 2005, MNRAS, 364, 1105 

Springel V., White S. D. M., Jenkins A., Frenk C. S., Yoshida N., 
Gao L., Navarro J., Thacker R., Croton D., Helly J., Peacock 
J. A., Cole S., Thomas P., Couchman H., Evrard A., Colberg J., 
Pearce F, 2005, Nature, 435, 629 
Springel V., Frenk C. S., White S. D. M., 2006, Nature, 440, 1 137 
Springel V., Wang J., Vogelsberger M., Ludlow A., Jenkins A., 
Helmi A., Navarro J. F, Frenk C. S., White S. D. M., 2008, MN- 
RAS, 391, 1685 

Stadel J., Potter D., Moore B., Diemand J., Madau P., Zemp M., 
Kuhlen M., QuiUs V., 2009, MNRAS, 398, L21 

Thomas P. A., Colberg J. M., Couchman H. M. P., Efstathiou G. P., 
Frenk C. S., Jenkins A. R., Nelson A. H., Hutchings R. M., Pea- 
cock J. A., Pearce F R., White S. D. M., 1998, MNRAS, 296, 
1061 

Thomas P. A., Muanwong O., Pearce F. R., Couchman H. M. P., 

Edge A. C, Jenkins A., Onuora L., 2001, MNRAS, 324, 450 
Tormen G., Bouchet F R., White S. D. M., 1997, MNRAS, 286, 

865 

van den Bosch F C, 2002, MNRAS, 331, 98 

Wechsler R. H., Bullock J. S., Primack J. R., Kravtsov A. V., 

Dekel A., 2002, ApJ, 568, 52 
White S. D. M., Rees M. J., 1978, MNRAS, 183, 341 
White S. D. M., Frenk C. S., 1991, ApJ, 379, 52 
White S. D. M., 1996, in Schaeffer R., Silk J., Spiro M., Zinn- 

Justin J., eds, Cosmology and Large Scale Structure Formation 

and Evolution of Galaxies, pp 349—1- 



Dynamical State of Dark Matter Haloes I 13 



